Pairing of Cooper Pairs in a Fully Frustrated Josephson Junction Chain 
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We study a one-dimensional Josephson junction chain embedded in a magnetic field. We show 
that when the magnetic flux per elementary loop equals half the superconducting flux quantum 
<f>o — h/2e, a local Z2 symmetry arises. This symmetry is responsible for a nematic Luttinger liquid 
state associated to bound states of Cooper pairs. We analyze the phase diagram and we discuss 
some experimental possibilities to observe this exotic phase. 
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During the last twenty years, Josephson junction ar- 
rays have proved to be very good tools to investigate 
classical and quantum phase transitions Recently, 
much attention has been payed to systems which display 
highly degenerate classical ground states S due to the 
presence of Aharonov-Bohm cages ||. Interestingly, a 
glassy vortex phase without disorder has been predicted 
for such two-dimensional (2D) structures H in agreement 
with experimental observations j5|. In the present work, 
we investigate the influence of quantum fluctuations on 
such systems in a ID model introduced in jq]. In this re- 
markably simple example, the huge classical degeneracy 
is a direct consequence of a local H-2 symmetry which is 
unbroken in the presence of quantum fluctuations. We 
show that this may stabilize an unusual nematic Lut- 
tiger liquid (LL) phase in which charge 4e bound states 
of Cooper pairs are the elementary objects. 

We consider the chain of loops shown in Fig. ^ em- 
bedded in a uniform magnetic field. We denote by <p 
the magnetic flux per elementary plaquette and we set 
7 = 2it(j)/4>Q where (j)Q = h/2e is the superconducting 
flux quantum. Each site of this lattice is supposed to 
be occupied by a superconducting island. A convenient 
description of the low-energy Hilbert space of this sys- 
tem involves local boson operators a'^b'^c^ (a n ,b n ,c n ) 
that create (destroy) Cooper pairs on the three types of 
islands of the lattice, respectively represented by black, 
grey and white circles in Fig. [l]. The system is described 
by the following Josephson coupling Hamiltonian : 



H, 



\{K + Cn + K-1 +e I7 c„_i) + h. c. (1) 



We first focus on the special value 7 = ir (half a flux 
quantum per loop). As shown in p], this Hamiltonian 
has, in this case, a single particle spectrum composed of 
three highly degenerate flat bands £0 = 0, e± = ±2£j. 
The corresponding eigenstates can be chosen as strictly 
localized (cage states) around each fourfold coordinated 



site (see Fig. [I]). This leads naturally to the notion of 
Aharonov-Bohm cages discussed in M . 




FIG. 1. The chain of loops and the three (nonnormalized) 
cage eigenstates corresponding to eq (left) and e± (right). The 
dashed line symbolizes the hopping term —tj e~ n . 

Let us introduce the set of boson operators A^ a n (A a>n ) 
that create (destroy) one Cooper pair with energy e a 
(a = 0, ±) in a cage state localized around the nth four- 
fold site. These operators can be simply expressed as a 
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linear combination of the operators ajj, &L c^, b' n _ 
only, whose coefficients are given in Fig. |l| so that we get 



Hi 



A 



(2) 



In this present form, Hj clearly exhibits a local U(l) 
symmetry. We shall now study the effect of boson-boson 
interaction on this symmetry. Therefore, we consider a 
real valued function n 1— > s n , and we construct a unitary 
operator U s defined by : UsAa^U^ 1 — e" 1Sn A a , n which 
commutes with H j . Using the precise form of the cage 
states, we easily obtain : 

UsCilanU^ 1 = ala n (3) 
UsblbnU- 1 = cos 2 (A„)&t&„ + sin 2 (A n )cJ l c„ + z n (4) 
Usclcnl/- 1 = sin 2 (A„)6j l 6„ + cos 2 (A„)cJ l c„ - z n , (5) 
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where we have set z n = isin(2A n )(&T Cn — cJj6 n )/2 and 
A n = (s n +l — s n )/2. From these transformation laws, 
we can readily see that any interaction term involving 
the bilinear operators a^a™, b^bn + c^Cn preserves the 
local t/(l) symmetry. Physically, this symmetry implies 
that the total number of bosons in each cage is separately 
conserved and the system remains an insulator. However, 
this symmetry is fragile since it is easy to find two-body 
interactions which break it. For instance, a Hubbard-like 
interaction term Sri( a n a «) 2 + (b n bn) 2 + (c}jC n ) 2 , has this 
effect which is manifested by the appearance of delocal- 
ized two- particle bound states discussed in Q . An excit- 
ing feature of this system is that this type of interaction 
still preserves a subgroup of the full £7(1) correspond- 
ing to a local Z2 symmetry. This subgroup corresponds 
to s n = [ir] for all n. With this restriction, it is easy 
to check that the operator b^bmb^bn + c^CmC^Cn com- 
mutes with U s for all (m, n). This local Z2 symmetry has 
an important physical consequence since it means that 
the parity of the total number of bosons in each cage is 
separately conserved. Therefore, if two-particle interac- 
tions lead to coherent transport through the chain for a 
many-boson system, quasi off-diagonal long range order 
may occur only for composite objects built with an even 
number of original bosons i. e., here, of Cooper pairs. In 
other words, a superconducting Josephson junction chain 
with this geometry and half a flux quantum per loop may 
realize a quasi-Bose condensate (in fact a LL) of charge 
4e composite bosons ! 

To discuss in more details the physics of this system, 
it is useful to rephrase these symmetry considerations in 
the language of quantum rotor mo dels. These offer the 
advantage of an intuitively simple classical limit defined 
from the phase of superconducting order parameter. For- 
mally, we introduce three phase fields 6 n ,ip ni \ n and their 
canonically conjugate fields He,n, n Vjn , n Xjn , which are 
related to the local Bose operators by : a\ L = e l9n , 

6jj = Tltp{n e ltPn , cjj = U]/,n e IX " . Assuming that the local 
particle number fluctuations are small, we get the quan- 
tum phase Hamiltonian |lj 

n 

Ej ^2 cos (°n ~ fn) + C0S(6>„ - \n) + 

n 

cos(6»„ - ) + cos(6>„ - Xn-l - 7) , (6) 

where Ec is the charging energy and Ej the Josephson 
coupling between islands. Note that the present mod- 
elling of capacitive effects is not meant to be very realis- 
tic, since for the sake of simplicity, we have not taken into 
account off-diagonal elements of the capacitance matrix. 
This choice corresponds to a local Hubbard-like interac- 
tion term between Cooper pairs. For convenience, we set 

V EqEj = 1. 



The classical ground state of H is easily obtained for 
any 7. Indeed, if we set x n = 9 n+ \ ~6 n — 7/2, and elimi- 
nate ip n and Xn, minimizing H is equivalent to minimize 
F(x n ) = - |cos (x n /2 + 7/4)| - |cos (x n /2 - 7/4) | for all 
n. As shown in Fig. |2| F has two local minima in x n = 
and x n = ir. For < 7 < n, one has F(0) < F(tt) and the 
classical ground state is unique (up to a global U(l) de- 
generacy). By contrast, for 7 = tt, one has F(0) = F(tt) 
so that, for a given plaquette, we obtain two degenerate 
ground states (up to a global translation of the phase 
variables) which are illustrated in Fig. S. 
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FIG. 2. Behaviour of F(x n ) for 7 = (solid line), 7 = 3vr/4 
(dashed line) and 7 = -k (dotted line). 

These states only differ in the sign of the superconduct- 
ing currents which circulate around the plaquette. For 
a chain made up of N loops, we thus get 2^ degenerate 
classical ground states up to a global translation of the 
phase variables. 




FIG. 3. Two possible classical ground states of H with dif- 
ferent chirality. 

This huge degeneracy is a direct consequence of a lo- 
cal Z2 symmetry of H . Note that H is not invariant 
under the full local U{1) group related to the Aharonov- 
Bohm cages. This occurs since the Josephson term in 
H may be written as a strongly non linear expression of 
the basic local Bose operators. For the Z2 transforma- 
tions, it is an easy task to translate the U s operators in 
the language of phase variables. An interesting local Z2 
transformation is provided by a kink in the s m 's. Let 
U n be the transformation defined by s m — for m < n 
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and s m = 7r for to > n. This transformation does not 
modify the phase variables for to < n whereas it shifts 
them by 7r if to > n and it permutes ip n and %„. Thus, 
its main physical effect is to change the currents flow- 
ing around the plaquette located between n and n + 1 
into their opposite value. From this description, we de- 
duce that starting from a given classical ground state, 
we may generate any other ground state by applying a 
finite sequence of such U n operators. We also see that, in 
the classical limit considered here, the local ~2-2 symme- 
try is spontaneously broken, yielding ground states with 
well-defined local circulating supercurrents. Note that 
U n also leaves most conjugate variables unchanged ex- 
cept H(p >n an d n Xj „ which are exchanged. As a result, 
we may add to H any term involving these conjugate 
fields without breaking the local ILi symmetry, provided 
the spatial symmetry between the tp n and Xn degrees of 
freedom is respected. Experimentally, this would require 
a tight control of offset charges since these may seriously 
alter an otherwise excellent geometrical symmetry of the 
chain. 

For real systems, it may become important to take into 
account quantum fluctuations of the phase variables, es- 
pecially when the superconducting islands are so small 
that their charging energy Ec can no longer be neglected 
in comparison to the Josephson coupling energy Ej . For 
a single loop and 7 = it, these quantum fluctuations have 
been shown, theoretically [Q and experimentally 
to induce tunneling between the two degenerate classical 
ground states shown in Fig. |^. The true quantum me- 
chanical ground state is therefore a macroscopic linear 
superposition of these two classical states and provides 
a simple example of a "Schrodinger cat" . For a system 
with N loops, eigenstates are classified according to the 
various irreducible representations of the local IL2 group 
which mix all the 2 N classical ground states. One of our 
next goals is to describe how quantum fluctuations lift 
the degeneracy among these representations, which is an 
artifact of the classical limit. 

At small g — yj Ec / Ej , the properties of the sys- 
tem are actually very similar to those of a quantum XY 
model. For small 7, we thus expect the infinite chain to 
be in a LL phase for g < g* (7) and in a gapful insulating 
(I) phase for g > g* (7). The transition at 0^(7) is of 
Berezinskii-Kosterlitz-Thouless (BKT) type fu|. Sim- 
ple spin wave calculations using the harmonic approx- 
imation of H around its classical ground state predict 
g* (7) = Vcos( 7 /4) 5 *(0) with fir* (0) = 7rV3/2. The main 
effect of the magnetic field, in this simple approximation, 



energy physics of this system. Therefore, instead of H, 
we now consider the following Hamiltonian : 



is thus to replace g by an effective g e ff — g / y/ cos(~f / 4) 
which controls all the correlation function exponents in 
the LL phase. 

To analyze the effect of the quantum fluctuations in 
the vicinity of 7 = tt where the additional local Z2 sym- 
metry emerges, it is convenient to eliminate the twofold 
coordinated islands to get a simple description of the low- 



^ 2 

n 

ecos(6>„ - 9 n+1 - 7/2)} 



-{cos[p(#„ - e n+1 - 7/2)] 

9 



(7) 



The parameter g' is provided by fitting the exponent of 
the correlation function (e 1 ( e " l ~ 9 ™- ) ) in the semi-classical 
regime with its value obtained with H in the harmonic 
approximation. This choice leads to g' = 2 5 / 4 3~ 1 / 2 g. 
The parameter e = 4|-7 — 7r| is determined from the energy 
splitting between the two local minima of the single loop 
potential energy. Finally, in our case, we have p = 2 
but we discuss thereafter the properties of Hxy for an 
arbitrary p. 

The Hamiltonian Hxy has a local TL p symmetry at 
e = 0, corresponding to the local transformations T a : 
®j l— * @j + 2iraj/p where a,j is an integer. The irre- 
ducible representations of this group are easily obtained 
in a basis which diagonalizes simultaneously the Ilea's. 
For a state \ip) such that U$ t j\ip) — lj\ip), where lj 
is an integer, we have T a \^p) — exp(iJ3j lj a,j)\ip). 
Writing lj = rrij + p nj with nij and rij integers and 
< nij < p — 1, we find that the set of m 3 's com- 
pletely specifies the irreducible representation of the lo- 
cal 2 p group. For each such representation, the action 
of the corresponding projector on the approximate gaus- 
sian ground state of Hxy produces a natural trial wave- 
function at least when j' < 1. We have computed the 
expectation value of Hxy on these states. Doing so, we 
noticed that 27r/p-tunnel processes occuring on different 
lattice sites are mostly uncorrelated. Neglecting com- 
pletely these correlations we get : 



(Hxy) 



g'L 2 _ Ce~f 
~2N~ V 



E 



cos 



2tt 



(to, - L/N) 



(8) 



up to a constant energy independent of the representation 
and to factors of order e~ 2 ? . In (js|), L = J2j h is the total 
angular momentum, C is a number close to 2tt 2 — 8 at 
small g' , and / ~ An/pg'. The ground state is therefore 
obtained by choosing the identity representation of the 
local Z p group (rrij = 0). 

Next, we see that the term proportionnal to e cou- 
ples different irreducible representations of the local Z p 
group. When p = 2 the action of this perturbation on the 
2 N low-energy trial states just discussed, is well described 
by a quantum Ising model in a transverse magnetic field : 



y \ n n / 



(9) 



where D is close to 1 for small g' . In terms of these Ising 
variables, the local 7Li symmetry corresponds to inter- 
changing the |+) and |— ) states on any given subset of 
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sites. It is therefore implemented by the a* operators. 
This model has a continuous phase transition (in the uni- 
versality class of the 2D Ising model) at its self-dual point 
which corresponds to e = (C/D) e~ 27r ^ g . Furthermore, 
it is easy to see that Hxy exhibits a BKT transition at 
7 = 7r (so e = 0) corresponding to the loss of quasi long 
range order for the nematic order parameter e 2W . In the 
harmonic approximation, this occurs at g(ir) — g*(n)/i. 

Gathering these informations obtained in various lim- 
its, we get the phase diagram drawn in Fig. ||. Besides 
the two familiar phases, namely the (I) phase at large 
g and the LL phase characterized by an algebraic order 
of the e lS order parameter, the most interesting result 
is the presence of a remarkable nematic Luttinger liquid 
(NLL) phase which may be viewed as a quasi-ordered 
condensate of pairs of Cooper pairs associated to the or- 
der parameter e . 
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FIG. 4 

g = \JE C /E 3 , and 7 
liquid. 



Schematic phase diagram as a function of 
2n(j}/(j)o- LL stands for Luttinger 



The dotted line of 2D Ising-type separates the two 
Luttinger phases. The Ising order parameter is vanish- 
ing in the NLL phase and builds up in the conventional 
LL phase where its square is proportional to the alge- 
braically decaying part of the e autocorrelation func- 
tion. The physical picture presented here has already 
been uncovered by Lee and Grinstein fll|| in the frame- 
work of a 2D classical XY model with squared cosine 
interaction. In their language, the BKT transition from 
the 2e LL phase to the I phase involves 2tt vortices of 
the 9 field, by contrast to it vortices between the NLL 
phase and the I phase. The nature of the multicritical 
point Q remains mysterious to us. Note that the cou- 
pling between XY and Ising-like degrees of freedom has 
been extensively studied |L2]-[l7|], already in the context 
of fully frustrated Josephson juntion arrays. Neverthe- 
less, in these works, the Ising order parameter has a very 
different nature since it describes the possible melting of 
the vortex lattice. It implies that the XY transition has 
to appear at lower temperatures than the Ising transition 



]1 2|— pL 4] . In our case and in the vicinity of 7 = it, the XY 
transition occurs for larger g than the Ising one. So, the 
effect of the Ising domain walls on the phase stiffness is 
very different in the two situations. 

To conclude, we have established the possibility of 
macroscopic condensation (in the sense of a LL) of charge 
4e objects in a Josephson junction chain for 7 ~ it. Ex- 
perimentally, it may be possible to detect this binding 
of Cooper pairs by connecting such a network to super- 
conducting leads. Indeed, we expect here a phenomenon 
analogous to Andreev's reflection : an ordinary Cooper 
pair of charge 2e entering the chain at low energy (com- 
pared to the gap between different representations of the 
local IL2 group ) will leave behind a pair of charge — 2e so 
that a charge 4e composite object may propagate along 
the chain. Another possibility is to close the chain into a 
large ring. In this geometry, we expect quantum oscilla- 
tions of the global current with respect to the magnetic 
flux across the ring with an elementary period </>o/2 as 
long as 7 ~ it. 
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